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Abstract 

Silent universes are studied using a "3+1" decomposition of the field equa- 
tions in order to make progress in proving a recent conjecture that the only silent 
universes of Petrov type I are spatially homogeneous Bianchi I models. The infi- 
nite set of constraints are written in a geometrically clear form as an infinite set 
of Codacci tensors on the initial hypersurface. In particular, we show that the 
initial data set for silent universes is "non-contorted" and therefore (Beig and Sz- 
abados [P) isometrically embeddable in a conformally flat spacetime. We prove, 
by making use of algebraic computing programs, that the conjecture holds in the 
simpler case when the spacetime is vacuum. This result points to confirming the 
validity of the conjecture in the general case. Moreover, it provides an invariant 
characterization of the Kasner metric directly in terms of the Weyl tensor. A 
physical interpretation of this uniqueness result is briefly discussed. 
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1 Introduction 



In the last few years there has been intense activity in the study of silent universes, 
which are dust spacetimes such that the fluid velocity vector is irrotational and the 
magnetic part of the Weyl tensor with respect to this vector vanishes. The concept was 
first put forward by S.Matarrese, O.Pantano and D.Saez in an astrophysical and 
cosmological context in order to describe structure formation in the universe. More- 
over, silent universes were interesting also from an exact solutions point of view because 
this class contains the Szekeres family of spacetimes ||^, which has been extensively 
analyzed in the literature (see |Q for a comprehensive review). The Szekeres class 
exhausts the silent universes of Petrov type D. Since silent universes are necessarily of 
Petrov type 0, D or I (Barnes and Rowlingson |Q) and the conformally flat case is well 
understood, the interesting new physics was in the Petrov type I subclass and a num- 
ber of relevant new results were expected. This fact triggered active research on silent 
universes @-in|. In particular, it was soon realized that the Einstein equations and 
the Bianchi identities, when written in a suitable orthonormal tetrad, decouple into a 
set of ordinary differential equations (which describe the time evolution of the space- 
time) and a set of constraint equations (which only involve derivatives along spatial 
directions). Most of the studies focussed in analyzing the evolution equations, mainly 
by making use of a dynamical system formulation of the Einstein field equations (see 
fl^ for a detailed account of the method). The existence problem of silent universes 
of Petrov type I was apparently settled in |TD[ where it was claimed that solving the 
constraint equations on an initial hypersurface is necessary and sufficient to obtain a 
silent universe (or, equivalently, that the time evolution of the constraints is automati- 
cally satisfied if the constraints hold initially). Unfortunately, this claim turned out to 
be untrue [l^ . So, the problem of how large the class of Petrov type I silent universes 
is had to be reanalyzed. In two independent works |]I^-|T5[ this issue was addressed. 
The analysis of the time-evolution of the constraints was performed in [|14| using a co- 
ordinate approach and in [|r^ using a tetrad approach. In both cases it was found that 
the successive time derivatives of the constraints become larger and larger expressions 
which are algebraically independent from each other. This led these authors to conjec- 
ture that silent universes of Petrov type I are extremely scarce and that the whole class 
reduces to some Bianchi models. Proving this conjecture turns out to be a very difficult 
problem because the successive constraints become huge just after a few time deriva- 
tives and they become unmanageable even for algebraic computing programs. Our aim 
in this paper is twofold. First, we perform a study of silent universes using the ADM 
splitting [|16] with respect to the hypersurface fl orthogonal to the fluid velocity. The 
main advantage of this method is that the full set of constraints can be written in a 
geometrically neat way as an infinite sequence of symmetric tensors satisfying the Co- 
dacci equation (i.e. a sequence of so-called Codacci tensors). Moreover, we show that 
the initial data set {Q, hat, Kab) satisfies the so-called non- contorted condition. R.Beig 
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and L.B.Szabados have recently proven that, locally, non-contorted initial data sets 
can be isometrically embedded in a conformally flat spacetime. Thus, each hypersur- 
face orthogonal to the fluid velocity vector in a silent universe can be locally viewed as 
a hypersurface in a conformally flat spacetime. We have not been able to exploit this 
fact fully yet, but we believe that an appropriate use of this result could prove essential 
to settle the question on the non-existence of spatially inhomogeneous silent universes 
of Petrov type I. However, this question is not analyzed here any further. Rather, we 
concentrate in the second objective of the paper. The great complexity of analyzing the 
set of constraints in the general case simplifles in the vacuum case (essentially because 
there are flve unknown functions instead of six). This allows us to prove, by making use 
of algebraic computing, that a vacuum silent universe of Petrov type I must be locally 
isometric to the Kasner spacetime (which is a one-parametric family of spatially 
homogeneous Bianchi I vacuum spacetimes) thus showing that the conjecture holds in 
vacuum. Although the concept of silent universe was originally put forward for dust 
spacetimes, the deflnition makes obviously sense for vacuum spacetimes as well. The 
only difference is that in vacuum there is no privileged timelike congruence, unlike in 
the dust case, and the deflnition must be modifled accordingly (see Deflnition |l|). Since 
Petrov type D vacuum silent universes are just the vacuum subclass of the Szekeres 
family, and vacuum conformally flat silent universes are, of course, locally Minkowski, 
our result provides a complete classiflcation of vacuum silent universes. This result is 
interesting at least in three respects. First, it suggests the validity of the conjecture 
in the general case. Furthermore, it provides an invariant characterization of the Kas- 
ner spacetimes without involving isometrics (the Kasner family is well-known to be the 
most general spatially homogeneous Bianchi type I vacuum metric). Finally, this result 
can also be interpreted physically as follows. The uniqueness result we prove states 
that non-trivial vacuum silent spacetimes must be of Petrov type D (or 0) (we consider 
Kasner as trivial due to its high degree of symmetry). Gravitational flelds of Petrov 
type D have often been considered analogous to Coulombian electromagnetic flelds in 
flat space. This analogy is, obviously, loose because the Petrov type of a spacetime is 
a local property, while characterizing an electromagnetic fleld as Coulombian involves 
conditions on the decay rate near spatial inflnity, which is a global property. Never- 
theless, accepting this analogy and considering Petrov type D gravitational flelds as 
"locally Coulombian", our uniqueness result could be interpreted physically as estab- 
lishing that vacuum gravitational flelds which are purely electric with respect to an 
inertial observer must be of Coulombian type (except for the Kasner metric, which is a 
special case). A similar result also holds in electromagnetism in flat space, where it is 
true that electromagnetic flelds which are purely electric with respect to an inertial ob- 
server must be Coulombian. However, this analogy should be taken cautiously because 
the terms "inertial" and "Coulombian" have a local meaning in gravity (corresponding 
to " geodesic irrotational observer" and "Petrov type D" , respectively) while they are 
global in flat space (meaning "orthogonal to hyperplanes" and "with decay rate as 
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r^^ at infinity"). Nevertheless, we believe that this analogy provides some physical 
explanation for the uniqueness result we have obtained. 

The paper is organized as follows. In section 2 we describe briefly the ADM formal- 
ism for an arbitrary energy-momentum tensor. This fixes our notation and conventions. 
In section 3 we adapt these equations to the silent universe case. In particular, we show 
that the evolution equations can be written as a system of ordinary differential equa- 
tions when suitable variables are chosen. Then, we concentrate in the study of the 
constraint equations. We prove that the whole set of constraints take the form of an 
infinite sequence of tensors {Sn)b which satisfy the so-called Codacci equation. The 
tensors depend algebraically on the second fundamental form Kab and the Ricci 

tensor of hat- Then, we show that the Codacci equation for (S'o)ab and {Si)ab imply 
that the initial data set is non-contorted and quote the result by Beig and Szabados 
mentioned above. In Section 4 we study the Codacci equation by rewriting it as a 
Pfaffian system for one-forms. This allows us to study the integrability conditions of 
this equation easily and show that {Sn)ab are symmetric and commute with each other. 
Using these results, the evolution equations are rewritten in terms of the eigenvalues 
of {Sq)1; and (^i)^. In section 5 we restrict ourselves to the vacuum case and we prove 
our main theorem, namely that an arbitrary vacuum silent universe of Petrov type I 
is locally isometric to Kasner. The computations necessary to show this result have 
been performed using Reduce. However, the proof is not straightforward and needs 
some discussion. Hence we have included the source code of the program, as well as 
the necessary explanations in Appendix A. Reading this Appendix in detail requires 
some knowledge of Reduce but we believe that the general idea of the proof can be 
understood also by the non-expert just by following the comments we have included. 



2 Preliminaries and basic results 



In this section we review briefly the ADM formalism |]16l and write down its basic 
equations. This will fix our notation and conventions. 

The study of silent universes was motivated, among other reasons, by the fact 
that the Einstein field equations split into a set of ordinary differential equations (the 
evolution equations) and a set of constraint equations. In order to see why this happens, 
let us start by writing down the "3-1-1" equations for an arbitrary energy-momentum 
tensor. Consider a smooth irrotational timelike congruence of curves on a spacetime 
(A^,(7)Qwith unit tangent vector u. Then, there exists an embedded hypersurface fl 
which is orthogonal to u, i.e. is a three-dimensional manifold and ipo : Q Ai 
is an embedding such that (Po{u) = (it is obtained by lowering the index to u and 

-'^For simplicity we assume the spacetime to be C°° although lower differentiability would suffice. 
The signature of the metric g is (— f , 1, 1, 1), the Levi-Civita covariant derivative is denoted by V and 
our sign conventions of the Ricmann and Ricci tensor follow jisl 
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the -k denotes the pullback, as usual). Let us choose an arbitrary point p & Q and 
an open, connected neighborhood Up C Q of p with compact closure. Let t denote 
the proper-time parameter of the congruence satisfying t = on ^o{Up). Then, there 
exists a positive real number to such that the integral curves of u crossing ipoiUp) exist 
for t E lo = {—to, to) C M and do not intersect (po{Up) again (also for \t\ < to). Let 
$j denote the one-parameter local group of transformations generated by u. The map 
ipt = ^t°'^o\up is an embedding ipt : Up ^ A4 and (pt{Up) is a spacelike hypersurface in 
M.. The set Mlp = [j\t\<to ^t{Up) is diffeomorphic to Up x /q and {J^p, g\Mp) is globally 
hyperbolic. Given an r-index covariant tensor field T on Aip we define a one-parameter 
family of covariant tensor fields on Up by Tj^...j^{t) = ip^ {T)j_^ j^, t G /q (tensors on 
Up carry Latin indices and tensors on Ai carry Greek indices). The definition of Lie 
derivative yields ip* (£^T)^.^...^-^ = dtTj^...j^.{t). Define as usual the acceleration and the 
deformation tensor by = u^V yU^ and T^^y = where = 6^ + u^u^ 

is the projector. Up can be endowed with a family of symmetric tensors hab{t) and 
Kab{t), t e Jo by 

habit) =ipti9)ab^ Kab{t)=vt{^)ab- 

For fixed t, hab{t) is a positive definite metric on Up. We denote the corresponding 
Levi-Civita covariant derivative by and the Riemann and Ricci tensors by Rabcd{t) 
and Rab{t). For objects at t = we will drop the argument, i.e. /lafe(O), Kab{0) and 
will be written simply as hat, K^b and D respectively (and similarly for other tensors 
on Up). The Gauss and Codacci identities relating the Riemann tensors of {M.p, qImp) 
and {Up, habit)) are Il9| 

^lin)abcd = Rabcdit) + KUt)KUt) - Kad{t)KUt), (1) 

^tinu))abc = DlKM-DlKacit), (2) 

where IZ is the Riemann tensor of {M.p, g\Mp) with all the indices lowered, and 7i{u) is 
the tensor TZin^a/s^y = u^Tisais-f Since the Ricci and the Riemann tensors are equivalent 
in three dimensions, we can contract (|l|) with h"''^{t) without loss of information. The 
result is best written by introducing the electric and magnetic parts of the Weyl tensor 
in iMp,g\Mp) as Eap = u^u'^C^.aup, = luf'u''r]^aTaC'"' ^p, ivap-ys is the volume 
form of the spacetime). Using now the Einstein field equations for an arbitrary energy- 
momentum tensor, we can decompose the Ricci tensor of {J^p, g\Mp) as 

TZafB = ^af3 + {P + P) Ua^/S + ^ (P " P) 9af3 - U^Qp - Upq^. 

where Yiap is symmetric and trace-free and IIq/j, are orthogonal to u. Let us define 
Eabit), Habit), Habit), Qait) , ^^(t), p(t) and p(t) as the pullbacks with respect to (pt of 
Ea/s, Hai3, Ilai3, Qa, o,j3, p and p respectively. The contracted Gauss identity reads 

Eabit) + \Tlabit) + \pit)habit) = Rabit) + Kit)Kabit) - Kacit)K'bit) = (3) 

= Sabit)+Sit)habit), 
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where K{t) = h''\t)Kah{t), S{t) = h''\t)Sab{t), Latin indices are raised with hab{t) 
and the symmetric tensor Sab(t) is defined through the second identity. The Codacci 
identity (0) can be written as 

V'dbc Hi{t) + \ ( h^t) q,{t) - habit) qrXt)) = DlKUt) - DlK,a{t) = Z^U^)- (4) 



where Zchaif) is defined by the second identity and r^*^^ stands for the volume form of 
{Up, habit)). Applying ip* to the trivial identity Cuga/s = ^ aUp + V/jMq, yields 

^ = 2KJf). (5) 

Finally, applying ip* to C^T.^^ = H^pLl^^ + a^a^ + h^h^V aap-h^uf^hJu^Rap^s (which 
is a well-known, and in any case easily verifiable, identity) we get, after raising one 
index and using @, 



dt 



-K^'.miit) + a\t)aXt) + Dy{t) - S',{t) + n\{t) - ^-^6^ (6) 



The geometrical identities are the so-called constraint equations and (^)-(||) are 

the evolution equations. 



3 ADM formalism for silent universes 

The evolution equations (H)-® are, in general, partial differential equations. The 
main property of silent universes is that, using an appropriate set of variables, the 
evolution equations become ordinary differential equations. To show this, we first 
write down the evolution equation for R%{t). Let r^^(t) denote the Christoffel symbols 
of habit). An easy consequence of (|) is dtT'^^{t) = DlKl{t) + ^/^(t). The identity 
dtR%c^{t) = Dl [dtT1^^{t)] — [dtT'l^{t)] becomes, after using the Ricci identity. 



dR%it) 

at 



D,Z''\ - DtZ''\ + KKl - KR% - ^RKl + 6% (^^RK - Tt{RK)^ 



(7) 



where TT{RK)(t) = R%{t)K\{t) and " |^" means that all the objects enclosed (including 
the covariant derivative) are to be taken at the value t. The crucial fact in (^ is that 
spatial derivatives act only on the tensor Zabc{t)- Hence, by imposing conditions on 
this object it is possible to obtain a system of ordinary differential equations for K%{t) 
and Rlit). The simplest possibility is to set Zabcit) = 0, which, from equation (^, 
is equivalent to qa(t) = and Habit) = 0. Moreover, (||) shows that some evolution 
equations for adt), p(t) and Ilab{t) are needed in order to obtain a closed system of 
differential equations. Again, the simplest possibility is that they vanish identically. 
This motivates the standard definition of silent universe. 
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Definition 1 A spacetime {A4,g) is called a silent universe if it admits a unit timelike 
vector field u with is irrotational and geodesic, the magnetic part of the Weyl tensor 
with respect to u vanishes and the Ricci tensor takes the form TZap = p (uaUf^ + ^ga/B^ , 
where p is an arbitrary smooth function (possibly zero). 

The condition that u is geodesic is superfluous in the dust case (p 7^ 0) because of the 
contracted Bianchi identities, but it must be required additionally in the vacuum case. 
The term "silent" stems precisely from the fact that the evolution equations become 
ordinary differential equations and hence, no influence from neighbouring points (apart 
from the one encoded in the initial data) arises during evolution. 

From now on, the spacetime {Ai, g) will denote a silent universe. Then, the Codacci 
identity becomes 

Z^at) = DiKUt) - DlKlit) = (8) 
and the evolution equation for K%{t) reads 

dK\{t) 



dt 



-{K\Ki + S\\. (9) 



Regarding the remaining evolution equation, it turns out to be more convenient to use 
S^it) instead of R%(t). A straightforward, if somewhat long, calculation using (0) and 
(i yields 



dt 



2S1KI + K^^Sl - 2SKI - 2KSI + 5^ {SK - Tr(Sir)) 



(10) 



where Ti{SK){t) = S^{t)K\{t). In order to derive this equation the following algebraic 
identity (which is valid for any 3x3 matrix A) has been used 



Tr(A2) Tr2(A) ^ , ^ .^^^^3^ _ ^^^(^^^^^^2) ^ ^^^3^^) 



A'' = ii{A)A^ + ^^^ tLMllA + l'^- ' ^ ^ 



Thus, the evolution equations for silent universes are indeed very simple (the deflnition 
was designed for this purpose). In return, the set of constraints become highly non- 
trivial. Indeed, the set of constraints Zabdf) = provide, after time differentiation, new 
constraints which must be satisfled identically. Let us describe this in detail. For the 
initial data problem, a triple {Up, hat, Kab) satisfying the Codacci constraint D^^a.K'^c] — 
must be given. Then, the tensor Sab is deflned via Sab + Shab = Rab + KKab — KacK^. 
Next, K°jj{t) and Sl{t) are obtained as the unique solutions of the ordinary differential 
system (|)-(|ig) with initial data Kl{0) = K% and S%{Q) = 3% Afterwards habit) can 
be found as the unique solution of (the right hand-side is already known) satisfying 
hab = hab{0). It is easy to show that the tensor S%(t) constructed a posteriori from 
hab{t) and Kab{t) coincides with the solution of (|1^) we started from. It remains the 
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check as to whether the Codacci equation for K%{t) is fulfilled for all t. A necessary 
condition is that all the time derivatives of Zabc{t) vanish at t = (whether this is also 
sufficient would require extra, non-trivial work). In [1^| and the time derivatives 
were discussed in a coordinate and tetrad setting respectively, and they were found 
to become increasingly large and intractable (the explicit expressions for the first few 
constraints appear in |2^) but almost nothing about the geometric structure of the 
successive constraints was revealed. The ADM formalism turns out to be more useful 
for this purpose. To that aim, let us find an expression for the commutation between dt 
and D^. Let NKt) be any one-parameter family of tensors on Up. A simple calculation 
yields 



d_ 

dt 



'd_ 
dt 



Lemma 1 Let {Ai,g) be a silent universe and {Up,hab{t),Kab{t)) constructed as in 
Sect. 2. Let A\{t) be a family of tensors on Up satisfying D*^y4''^](t) = Vt. Then the 
family of tensors 

BMt) = d,A\it) + K',{t)Ai{t) 
also satisfies D\^B\p) = Vt. 

Proof. Taking the time derivative of A^^j (t) = and using the commutation formula 
dnl) we find 



d 

0= —[DwA 
dt 



where Dj^^K^^-^it) = was used in the second equality and Dj^^A^^-^{t) = in the third 
one. □ 
This lemma, together with the constraint equations (^ suggests defining the fol- 
lowing sequence of one-parameter families of tensors on Up 

(SoUt) = -KUt), (SMt) = I [(Sn-im)] + K\{t){Sn^,fXt\ nen. (12) 

Lemma |I] shows that {Sn)\{t) satisfies D\^{Sn)\^{t) = 0, Vn G NU {0},Vt. Hence, the 
successive time derivatives of Zahcii) = at t = can be written as (S'^)'^^] = 0, 
Vn G NU {0} and these are the full set of constraints that the initial data {Up, hat, Kah) 
must satisfy. 

Notice that the evolution equation (|^) implies {Si)\{t) = S%(t). Furthermore, if 
follows from (^ and ([T0|) that each (5'„)"^(t) depends algebraically on K%{t) and S%{t). 
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Moreover, it follows easily that this dependence is polynomial and that if we decide 
that K%{t) carries a degree equal to one and Slit) carries a degree equal to two, then 
{Sn)\{t) is homogeneous of degree n + 1 (i.e. it consists of a sum of terms of degree 
n + 1). Unfortunately, obtaining an explicit general formula for {Sn)\{t) in terms of 
K1{t) and S1{t) seems to be a difficult task. The first two constraint equations read 
simply 

In particular, this shows that the initial data set {Up, hab, Kab) is non-contorted (see 
for a definition). In that paper, the authors prove that any non-contorted initial data 
set can, locally, be isometrically embedded in a conformally fiat spacetime. Hence, each 
hypersurface orthogonal to the fiuid fiow in a silent universe can be locally embedded 
in a conformally fiat spacetime. This already gives substantial geometric information 
about the initial data set (and only the first two constraints have been used!). We 
believe that this fact can provide the key for proving the conjecture that silent universes 
of Petrov type I must be spatially homogeneous (in terms of the initial data set, being 
spatially homogeneous is equivalent to {Up, hab) being fiat and Kab being covariantly 
constant). However we have not been able to exploit this fact fully yet and the question 
remains under investigation. 

Summarizing, in this section we have obtained the evolution equations for the silent 
universes as ordinary differential equations in terms of the tensors K\{t) and S%{t) 
and we have written the full set of constraints in the form D*^ (S'n)*c](^) ~ ^"^^ ^ 
collection of tensors {SnYb{t) which depend algebraically on K%{t) and S%{t). Since all 
these tensors satisfy the same equation, the next section is devoted to study it in some 
detail. 



4 Consequences of the Codacci equation 

Let us start with a standard definition 

Definition 2 Let (V, 7) he a (pseudo-)Riemannian manifold and D the Levi-Civita 
covariant derivative. A symmetric tensor field Qab is called a Codacci tensor iff it 
satisfies the Codacci equation 

DlaQ%] = 0. (13) 

Codacci tensors have received considerable attention in the mathematics literature 



see 



pm for an account) mainly because of the crucial role they play in isometric 
embeddings of Riemannian manifolds into Euclidean manifolds. In our context they 
are interesting because the full set of constraints for silent universes take the form of an 
infinite set of Codacci tensors (5'„)°j in the initial data set {Up, hab, Kab) (the symmetry 
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of {Sn)ab will t)6 shown later). Let us now consider the integrability conditions of the 
Codacci equation (|13]). To that aim, we will rewrite this equation as a differential 
system for one-forms (this way of writing the Codacci equation is new, as far as we 
know). Let Qab be a Codacci tensor on (V,7). Choose an arbitrary cobasis {0"'} on 
this Riemannian manifold and denote by ujI the connection one-forms in this basis. 
The torsion-free condition is dO"" + a;"^ A 9^ = 0. Denoting by the components of 
the Codacci tensor Q (with one index raised) in the basis {0"^} and its dual {ca}, we 
can define three one-forms by = Q1 6^ . A simple calculation gives 

dQ'' + AQ'= {d^cQ\^ e'Ae' = 0. (14) 

Thus, the Codacci equation for Qab is equivalent to dQ"" + lvI A = 0. To study 
the integrability conditions of this system we only need to take its exterior derivative, 
which gives 

n% AQ^ = o (15) 

where are the curvature two-forms n\ = duj'\ + a;" A oj'^f,. Let us now assume that 
V is three-dimensional (this is the relevant case for silent universes) and define the 
so-called Ricci one-forms by Pa = Rab {Rab are the components of the Ricci tensor 
in the triad {0°"}). The following identity holds in three dimensions 

Q"^ = P" A0b- PbAO" --RO" A 0b, 

where indices are raised and lowered with hab = h{ea, Cb) (this identity just states the 
well-known equivalence between the Riemann and Ricci tensors in three dimensions). 
The fact that Qab is symmetric can be written as A 6a = so that the integrability 
conditions ( |T5|) become Pb A = or, in index notation, 

R%Q',-Q%R', = 0. (16) 

This equation states that R^ and Q^f, commute when viewed as linear maps on the 
tangent space TgV, g G V. Let us now put (V,7) = (Up, hab)- Applying ([T6|) to the 
Codacci tensor we obtain, after using the definition of 5**^ (|^), 

SIK", - KIS\ = 0. 

Since all (SnYb polynomials in K^f^ and S% we obtain {SnYbi^"i)\ — {SmYbi^n)\ = 0, 
yn,m G N U {0}. In particular, each tensor {Sn)l commutes with K^. The defining 
recursion (0) then implies that {Sn)ab Vn G N are symmetric and hence Codacci 
tensors. 

Obviously, the same arguments hold for an arbitrary value of t G Iq. Thus {Sn)%{t) 
are Codacci tensors in {Up,hab{t)) and they commute with each other (for fixed t). 
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Therefore, there exists an orthonormal basis of vectors {ea{t)} in {Up, hab{t)) such that 
{SnYbif), 'in diagonahze simultaneously. Denoting by {0 {t)} the dual cobasis, we have 

Kit) = Xlit) e\t) ® ee(t), Sit) = Xlit) e\t) ® e,{t), (17) 

c=l c=l 

where Kit), Sit) denote the tensors with indices KUt) and S^t), and Ao(t) and \\it) 
are their respective eigenvalues. Writing equation ( ]T0| ) in this frame gives 

d\1 




dt 



+ "iX^Xl - 2S\l - 2K\l + SK~ TiiSK) 



(1^ 



where the Einstein summation convention has been suspended. Until here no restriction 
on the Petrov type of the silent universe has been imposed. However, silent universes 
of Petrov type D and are well-understood and attention can be restricted to Petrov 
type I. Of course, the Petrov type of a spacetime need not remain constant everywhere 
but it is well-known that the set of points with Petrov type I is open (see e.g. p2| ) 
and therefore a submanifold, so we can assume Petrov type I everywhere without loss 
of generality (because only local properties are considered in this paper). For silent 
universes, Petrov type I is equivalent (see e.g. p) to Eabit) having three different 
eigenvalues (of course only two of them are linearly independent because of the trecefree 
condition of Eabit)). Eq. (|^) shows that this is also equivalent to Sabit) having three 
different eigenvalues. Thus, (|T8|) for a 7^ 6 gives 



which allows us to rewrite equation (^ as 



(Xlit) + [X^') . (20) 



dt 

We now quote without proof a lemma due to Barnes and Rowlingson |^ which we will 
use in the vacuum case. The proof can be easily rewritten in our formalism by using 



the Codacci equation in the form (14) and the relations (p!8D ,(p!9D and (|20|). 

Lemma 2 (Barnes and Rowlingson j^) Let iAi,g) be a silent universe of Petrov type 
I and iUp, habit), Kabit)) constructed as in Sect. 2. Fix an arbitrary point q G Up. 
Then, 

1. the subset of Iq where K\qit) has three different eigenvalues is dense in Iq. 

2. Fix an arbitrary value t E Iq. Each one-form in the cobasis {0 (t)} introduced 
above is integrable. 

3. There exists a coordinate system {x,y,z} in Up in which Kit) and Sit) are 
diagonal and such that the metric habit) takes the form ds^it) = Ait, x, y, z)dx^ + 
Bit, X, y, z)dy'^ + C(t, x, y, z)dz^ . 
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5 Vacuum silent universes of Petrov type I 



Vacuum silent universes are characterized by p = which, from equation (|^), is equiv- 
alent to S{t) = 0. Hence, the evolution equation (|TH|) takes the simpler form 



dXljt) _ 
dt ~ 

We need the following lemma 



3A^AS - 2\1K - Ti{SK) 



(21) 



Lemma 3 Let {Ai,g) be a vacuum silent universe of Petrov type I and construct 
{Up, hab{t) , Kabit)) as in Sect. 2. Fix an arbitrary point q G Up. Let f^q{t) denote one 
of the eigenvalues A^|g(t), A = 0,1 defined in (|7?[j. Then, the set {t G /q; Pg(^) 7^ 0} is 
dense in Iq. 

Proof. Let first ^q(t) be one of the eigenvalues of S. We can choose fJ^qit) = X\\q(t) 
without loss of generality. Suppose that the lemma does not hold, i.e. that fiqit) 
vanishes for t G /i where Ii is open and non-empty. Then, equation (^I]) implies 
Tr{SK)\q{t) = Vt G h, which, after using S{t) = 0, implies (Ag(t) - Xlit))\q = 
Vt G /i. This is impossible from the first conclusion of Lemma |^. Let now fiq{t) be one 
of the eigenvalues of K{t), say Ag|g(t). The claim of the lemma follows easily because 
if fig(t) vanished on a non-empty open set Ji, then ( pO| ) would imply Xl\q(t) = on Ii, 
which we have just shown to be impossible. □ 
This lemma combined with Lemma ^ shows that there exits an open dense subset 
W G Uq X Iq where the eigenvalues Ag(t) are non-zero and mutually distinct and that 
the same holds for Xi{t). Then, the following parametrization exists on W 

Xl{t)=w\t, Xlit) = w{l + u)\t, Xlit) = w{l + v)\t, 
X\{t) = w^s{k + l)\t, Xlit) =w^sik-l)\t, Xlit) = -2whk\t, 

where w, u, v, s and k are nowhere vanishing scalar functions on W . Furthermore, the 
combinations u — v, v+1, u + 1, S/c — 1, 3/c + l, k — 1 and k + 1 are also nowhere zero on 
W (these statements just translate the fact that Sit) and Kit) have three different and 
non- vanishing eigenvalues on W). This type of parametrization is convenient for the 
algebra computing calculations we describe in Appendix A. The evolution equations 
(pop and (|2T| ) take the following form on W 



dfW = —w"^ isk + s + 1) , dtu = w (usk + us + 2s — — , 
dtV = w ( ksv + 3ks -|- (f + 1) (s — v)^ , dfk = w ( —ku + 2kv u + -k'^u 



2 2 

u 3 

dts = ws[2sk-l + 2s 2v - -ku ) . (22) 



We are now in a position to prove our main theorem. 
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Theorem 1 Let {Ai,g) be a vacuum silent universe of Petrov type I. Then the space- 
time is locally isometric to a Kasner spacetime. 



Before proving this resuh, let us recaU that the the Kasner spacetime is defined as 
the manifold x endowed with the metric 

ds^ = -dt^ + t^P^dx^ + t^^-'dy"^ + t^P'^dz^, t > 0, -oo < x, ^ < +00 (23) 

where Pi,P2,P3 G K and satisfy the relations P1+P2+P3 = pI+pI+pI = 1- Hence each 
element of the family is parametrized by one real parameter. This spacetime contains 
a three-dimensional abelian isometry group acting transitively on the spacelike hyper- 
surfaces t = const. Is is thus a Bianchi type I vacuum spacetime. Furthermore, it is 
well-known (see e.g. [|l^) that any spacetime admitting an abelian three-dimensional 
Lie algebra of Killing vectors which span spacelike hypersurfaces must be locally iso- 
metric to Kasner (a global isometry requires further conditions on the topology of the 
spacetime). Hence, the Kasner family can be locally characterized by the existence 
of these Killing vectors. Theorem ^ provides another characterization for the Kasner 
family directly in terms of the Weyl tensor and without involving isometrics, namely, a 
vacuum spacetime which is of Petrov type I and silent must be locally isometric to Kas- 
ner. Obviously a global isometry cannot be expected because the Petrov type I and the 
silent conditions are purely local and they place no restriction on the global topology of 
the spacetime (which could be, for instance, E x T'^ , where is the three-dimensional 
torus, endowed with the metric (p^)). 

The proof of Theorem |l| involves combining several constraint equations in order to 
show that they lead to incompatibilities unless the spacetime is very simple, namely a 
homogeneous Bianchi model. Since the expressions involved soon become very large, 
the proof would be impossible (using this direct method) without employing algebraic 
computing. However, the proof can be followed in exactly the same way as an ordinary 
proof can. The only requirement is some knowledge of Reduce, which is the algebraic 
computing program we have used. Even for those who are not familiar with Reduce, 
we believe that the general idea of the proof can be followed with relatively little effort. 
Hence, the Reduce program and the necessary explanations on the logic of the proof 
have been included in Appendix A (the meaning of the commands is not explained, for 
the interested reader we recommend the introduction to Reduce by M.A.H.MacCallum 



and F.Wright fH). 

Proof of the Theorem. 

Let us consider a point p E A4 and construct {Up, habit), Kabit)), t E Iq a.s described 
in Sect 2. Using Lemma 0, the metric hab{t) in Up can be written as 

ds'^{t) = A(t, X, y, z)dx^ + B{t, x, y, z)dy'^ + C{t, x, y, z)dz'^. (24) 

In Appendix A we proof that, under the assumptions of the theorem, the functions 
A, B and C depend only on t. Then, the spacetime metric in {M.p, g\Mj,) can be 
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reconstructed from habit) in the usual way to give ds'^\Mp = —dt'^ + A(t)dx'^ + B(t)dy'^ + 
C{t)dz'^. Thus, the spacetime admits three commuting Kilhng vectors which span 
spacehke hypersurfaces. Hence, the metric must must take the form (p3D for some values 
of pi, p2 and ps. Conversely, a simple calculation shows that the Kasner spacetime is 
a silent universe of Petrov type I. This concludes the proof of the theorem. □ 
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Appendix 

This Appendix contains the bulk of the proof of Theorem |l|. The aim is to prove 
that the one-parametric family of metrics (|2^ ) on Up (for vacuum silent spacetimes of 
Petrov type I) can be written in the form ds^{t) = A{t)dx'^ + B{t)dy'^ + C{t)dz'^. In 
order to prove this, we use the algebraic computing program Reduce. In this Appendix 
we include the code lines, which are written using the sans serif font, and we explain 
the logic of the proof (this is written in plain text). A few key output results^, which 
are required to follow the proof, are also included. These will be written in italic 
characters. 

load-package groebnerS on gcd$ 

depend w,t,z$ depend u,t,z$ depend v,t,z$ depend s,t,z$ depend k,t,z$ 

Only the dependence on t and z of the functions is prescribed. Of course, all 
of these functions depend on y and x as well, but we will not invoke any equations 
containing partial derivatives with respect to x or y. This may seem inadequate because 
we are not exploiting all the available equations. However, the strategy of the proof 
relies on showing a certain property for the z variable (the vanishing of some associated 
Christoffel symbols) and then use the property that the z variable is on the same footing 
as the X and y variables (there is an obvious symmetry between x< — >y< — >z 
in the metric (0)) in order to imply the vanishing of some other Christoffel symbols. 
Hence concentrating only on the z variable is sufficient. 

Let us now introduce the evolution equations (|2^ ) and the set of Codacci tensors 

s(0,l):=-w$ s(0,2):=-w*(l+u)$ s(0,3):=-w*(l+v)$ s(l,l):=w^2*s*(l+k)$ 
s(l,2):=w^2*s*(k-l)$ s(l,3):=-2*w^2*s*k$ 

^In Reduce, a code line ending with $ does not produce any written output (of course the operation 
is performed). If the hne ends with ";" the result is written in the screen (or the standard output). 
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dwt:=-w''2*(s*k+s+l)$ dut:= w*(u*s*k+u*s+2*s-u''2-u)$ 
dvt:= w*(v*s*k+v*s+s+3*k*s-v"2-v)$ 
dst:= w*s*(2*s*k-l+2*s-u/2-2*v-3/2*k*u)$ 
dkt:= w*(-k*u+2*k*v-u/2+3/2*k"2*u)$ 

df(w,t):=dwt$ df(u,t):=dut$ df(v,t):=dvt$ df(s,t):=dst$ df(k,t):=dkt$ 
d:=6$ for l:=2:d do for i:=l:3 do s(l,i):=df(s(l-l,i),t)-s(l-l,i)*s(0,i)$ 

The symbol s{l, i), where / G N U {0} and i = 1,2, 3, is used to denote the i — th 
diagonal coefficient of the Codacci tensor {Si). Let us introduce the Codacci equations 
(we denote the Christoffel symbols F*^ of the metric (p^) by gam(i,j,k), k = 1,2, 3). 

for l:=0:d do << cdz(l,l):= df(s(l,l),z)+gam(l,l,3)*(s(l,l)-s(l,3))$ 

cdz(l,2):= df(s(l,2),z)+gam(2,2,3)*(s(l,2)-s(l,3)) >> $ 
df(w,z):=rhs first solve(cdz(0,l),df(w,z))$ df(u,z):=rhs first solve(cdz(0,2),df(u,z))$ 
df(k,z):=rhs first solve(cdz(l,l),df(k,z))$ df(s,z):=rhs first solve(cdz(l,2),df(s,z))$ 
df(v,z):=rhs first solve(cdz(2,l),df(v,z))$ 
for i:=0:(d-3) do r(i):=num cdz(i+3,2)/w'^(4+i)/s$ 
for i:=0:(d-4) do n(i):=resultant(r(0),r(i+l),gam(l,l,3))/2 $ 

Our next aim is to analyze the case in which gam(2,2,3)= on a non-empty open 
subset W <ZW and gam(l,l,3) is nowhere zero on W . 

gam(2,2,3):=0$ 

pro:=r(0)/gam(l,l,3)$ proz:=df(pro,z)/gam(l,l,3)$ prozz:=df(proz,z)/gam(l,l,3)$ 
nl:=-2*v"2+12*v*(k-l)-9*k"2+30*k-13$ n2:=9*k-4*v-7$ 
fa ctorize { prozz- n 1 * pro- n2 * proz} ; 

{18,-v,v,3k+ l,3k+ l,k - 1} 

But this is impossible because no element in this list can vanish anywhere on W. 
This concludes this case. Next, we consider the situation in which gam(2,2,3) is non- 
zero almost everywhere on W (studying this case constitutes the heart of the proof). 
Possibly after restricting ly to a smaller dense open subset, we can assume gam(2,2,3) 
7^ everywhere on W (from now on we will restrict W to an open dense subset thereof 
whenever necessary and without further notice). 

clear gam(2,2,3)$ poll:=-n(0)/gam(2,2,3)$ pol2:=n(l)/gam(2,2,3)+5*u*poll$ 
ll:=22*( 585*k"2*s - 192*k*s - 25*s+360*k*u)$ 
I2:=-I989*k^2*u - 2712*k*u - 7008*k - 275*u$ 
pol3:=-264*k*n(2)/gam(2,2,3)+ll*poll+l2*pol2$ 

A number of particular cases must be analyzed. Let us start by studying what 
happens when the polynomial Cq, defined immediately below, vanishes on a non-empty 
open subset Wi C W. 

Subcase 0: Cq = identically on Wi. 

cO:= - 27*k''4*u''2 - 72*k"3*v*u + 36*k''3*u"2 + 96*k''2*v"2 - 96*k''2*v*u 
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+42*k'^2*u'^2 - 24*k*v*u + 12*k*u'^2 + u'^2$ 
c0t:=df(c0,t)/2/w$ c01:=resultant(cOt,poll,s)$ 

rg:=resultant(c0,c01,v)/(3*k+l)'^12/(3*k-l)'^12/kn2/3131031158784/u'^14$ 

rg+(-117*k^4+30*k^2-l )^2; 



So, — llT/c"^ + 30A;^ — = 0. In particular k is constant and its derivative is 
identically zero. Then 

part(groebner({dkt/w*2,-117*k'^4+30*k'^2-l,c0}),4); 

which is impossible because u cannot vanish anywhere on 1^. So we can assume 
Co 7^ everywhere on W and divide out this factor whenever necessary. This finishes 
Subcase 0. We continue with the generic case. 

v:=u*l$ ul:=resultant(poll,pol2,s)/(3*k+l)/(3*k-l)/(u-v)/v/u'^7/c0/324$ 
u2:=resultant(poll,pol3,s)/c0/162/v/u'^6/(u-v)$ 

By construction, the polynomials u\ and U2 depend only on A;, I and u. 

ql:=sub(k=2,l=2,ul)$ length coeff(ql,u); length coeff(ul,u); 

5 5 

q2:=sub(k=2,l=2,u2)$ length coeff(q2,u); length coeff(u2,u); 

7 7 

The two polynomials qi and q2 depend only on u. The two last output results show 
that qi and Ui are of the same degree as polynomials in u and the same happens for 
q2 and U2. 

groebner{ql,q2}; 
{1} 

Hence q\ and q2 have no common zeros. Prom this we can conclude that the 
resultant of Ui and U2 with respect to u is not identically zero. Indeed, if it were zero, 
this would mean that the polynomials Ui and U2 have common solutions for u at any 
value of k and /. In principle, it could happen that the common solution for u tends to 
oo when we approach k = 2 and / = 2, but this can also be excluded because qi and ui 
are polynomials in u of the same degree (and similarly for q2 and 1*2). So, there should 
exist at least one finite value of u which solves qi — Q and g2 = simultaneously. But 
we have proven that qi and q2 have no common zeros. Thus, the resultant of Ui and 
U2 with respect to m is a non-identically vanishing polynomial in / and k. The reason 
why we do not prove this by direct calculation is that the polynomials are already very 
big and performing the full resultant takes far too long. 

So, k and I must belong to the set of zeros of a polynomial in k and Z, which 
we denote by a (we view cr as as a subset of the {/c, Z} plane). Next, we prove that 
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restricting the polynomial ui to any point on a provides a non-trivial polynomial in 
u (i.e. non- constant). Similarly, we show that P0I2 (which depends on k,l,u,s) when 
restricted to any point on cr is a non-trivial polynomial in s. These two facts imply 
that {k,l,u,s} must lie on a one-dimensional manifold. So, let us first prove that ui 
is a non-trivial polynomial in u at any point on a. 

jl:=lterm(ul,u)/u^4$ j2:=lterm(u2,u)/u^6$ 
j3:=sub(u=0,ul)/(3*k+l)/(3*k-l)$ 
j4:=sub(u=0,u2)/(3*k+l)/(3*k-l)/414720*u'^2/k'^2/c0$ 
g:=groebner{j3,j4}$ g2:=second(g)/r2/(l-l)^2$ 

yl:=resultant(jl,g2,l)/(k-l)^3/(k+l)^5/kn2/(3*k+l)^7/(3*k-l)^7/471859200$ 

y2:=resultant(j2,g2,l)/3774873600/(3*k-l)^6/(3*k+l)^6/(k+l)^5/(k-l)^3/kn4$ 

groebner{yl,y2}; 

U} 

which proves that ui is nowhere identically zero on a. Since a constant non-zero 
value is also impossible (recall that -ui = is a consequence of the Codacci equations) 
the claim follows. Regarding P0I2, it suffices to notice that the leading term (in the 
variable s) of this polynomial is 

lterm(pol2,s)$ on factor$ ws; off factorS 

-44{3k + lf{3k-lfs^ 

which is nowhere zero on W. Hence, u, k, v, s lie on a curve. Therefore, there exists 
a non-zero vector field a = aid^ + a2dt (oi and 02 unknown functions on W) which 
annihilates k,s,u,v. Consider first the possibility that ai vanishes on a non-empty 
open set, i.e. suppose that k, s, u, v do not depend on t on that open set. Then 

v:=rhs first solve(2*dkt/w,v)$ s:=rhs first solve(2*k*dst/w/s,s)$ 

third groebner(num dut/w,num dvt/w,num poll/u,num pol2,u,k)/(3*k-|-l)"2/ 

(3*k-l)^2/(k+l)^3/k; 

4k^ — 1 Hence k is constant. 

first groebner({num dut/w,num dvt/w,4*k"2-l},{u,k}); 

{8k + 3u + 8, ~ ^} So, It is also constant. In particular, its derivative with 
respect to z vanishes. 

first groebner{df(u,z),r(0),2*u-F8-F8*k,4*k'^2-l}; clear s$ clear v$ 

gam{2,2,3) 

This is impossible in the case we are analyzing. So, we can assume without loss 
of generality that the vector a takes the form a — — b ■ gam(2, 2, 3)/wdt for some 
unknown function b. This provides a new set of equations (we also rewrite gam(l,l,3) 
= /*gam(2,2,3) where / is an unknown function) 

clear v$ gam(l,l,3):=gam(2,2,3)*f$ 

pr(l):=2*(df(k,z)-b*gam(2,2,3)*df(k,t)/w)/gam(2,2,3)$ 

pr(2):=2*(df(s,z)-b*gam(2,2,3)*df(s,t)/w)/gam(2,2,3)/s$ 
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pr(3):=(df(u,z)-b*gam(2,2,3)*df(u,t)/w)/gam(2,2,3)$ 

pr(4):=(df(v,z)-b*gam(2,2,3)*df(v,t)/w)/gam(2,2,3)$ 

pr(5):=r(0)/gam(2,2,3)$ 

for i:=l:4 do x(i):=resultant(pr(l),pr(i+l),f)$ 

zl:=-resultant(x(l),x(2),b)/(k-l)/(3*k+l)/2$ 

z2:=(-2*resultant(x(l),x(3),b)/(k-l)/(3*k+l)/2+(3*k+l)*zl)/(3*k-l)$ 

z3:=-resultant(x(l),x(4),b)/(k-l)/(3*k+l)/2$ 

fl:= ( 9*k'"2*u - 24*k*u + 30*k*v - 9*u + 6*v)/2$ 

f2:= ( -9*k'"2*u - 6*k*u - 36*k*v - 36*k + 3*u + 12*v + 12)/2$ 

z3:=(z3-fl*zl-f2*z2)/(3*k-l)$ 

Before following with the general case, we must study a few particular cases. They 
are characterized by the vanishing of one of the following expressions on a non-empty 
open subset W2 CW 

cl:= 2*k*v + 3*k + 2*v + 1$ 

c2:= 2*k*v+3*k+l$ 

c3:= - 3*k'"2 + 2*k*v + 2*k + 2*v + 1$ 

c4:=- 4*k"3*v + 3*k"3 - 10*k"2*v - 17*k"2 - 4*k*v - 3*k + 2*v + 1$ 
Subcase 1: ci vanishes identically on W2- 

clt:=2*df(cl,t)/w$ v:=rhs first solve(cl,v)$ clt:=clt/(3*k+l)*(k+l)'^2$ 
u:=rhs first solve(clt,u)$ rl:=-z2*8*(k-l)'^2/(k+l)/(3*k+l)$ 
rl-(2*k*s+2*s-l)^2*(5*k-l)^2; $ 

So, either k — 1/5 or 2ks + 2s — 1 = on some non-empty open set. k — 1/5 
is impossible as the following calculation shows 

k:=l/5$ groebner{dkt/w,poll}; clear k$ 

{1} Thus, 

s:=l/2/(k+l)$ factorize first groebner{poll,pol2}; clear v$ clear s$ clear u$ 

{3k + l,3k + 1} 

This is also impossible and therefore Ci cannot vanish on W2- As usual, we can 
assume that it is nowhere zero on W. 

Subcase 2: C2 vanishes identically on W2. 

c2t:=df(c2,t)/w*2$ gr:==groebner({zl,c2,c2t},{s,u,v})$ 
s:=rhs first solve(first gr,s)$ u:=rhs first solve(second gr,u)$ 
v:=rhs first solve(third gr,v)$ 
length coeff(num poll,k); clear u$ clear v$ clear s$ 

{8} 

Thus, k must by constant and hence s, v and u are also constants. In particular they 
do not depend on t and we have seen above that this is impossible for gam(2,2,3)7^ 0. 
Subcase 3: C3 vanishes identically on 1^2- 
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c3t:=df(c3,t)/w$ gr:=groebner({c3,c3t},{s,v})$ 

v:=rhs first solve(second gr,v)$ s:=rhs first solve(first gr,s)$ 

cc:=factorize (num zl/(3*k+l))$ 

length cc; 

2 

second groebner({first cc,z2},{u,k})$ on factorS ws; off factorS 

{3k^ + l){k^ + l){3k + l){3k - l){k - 1) 

second groebner({second cc,poll,z2},{u,k})$ on factorS ws; off factorS 
{3k + l){3k- lY{k- 1) 
clear uS clear vS clear sS 

This gives a contradiction to C3 = on hence C3 is non-zero everywhere on W . 
Subcase 4: C4 vanishes identically on W2 

c4t:=2*df(c4,t)/w$ groebner(coeff(c4,v)); 
{i } This shows that we can solve for v 
v:=rhs first solve(c4,v)$ 
groebner({zl,z2,5*k"2-l,c4t}); 

{ i } We can solve s in c4t 

s:=rhs first solve(c4t,s)$ length coeff(resultant(nunri zl.num z2,u),k); 

52 Hence, k must be constant and its derivative must be identically zero. Then 

second groebner({dkt,zl,z2},{u,k})/num(v)S on factorS ws; off factorS 

clear v$ clear s$ 

{k^ + l){3k- 1) 

Again a contradiction. This concludes the last particular case. 
Generic Case: We can now consider the generic situation (we divide by Ci, C2, C3 
and C4 whenever necessary) 

h(l):=resultant(zl,z2,s)/cl/2$ h(2):=resultant(zl,z3,s)/cl/2$ 
h(3):=resultant(zl,poll,s)S 

t(l):=resultant(h(l),h(2),u)/9216/k'^4/(k+l)'^2/(3*k-l)'^2/(3*k+l)'^2/c2/c3'^2/ 
c4'^2/(k-l)/(v+l)/v$ 

t(2):=resultant(h(l),h(3),u)/442368/k'^6/(3*k-l)'^4/v'^3/(3*k+l)'^3/c3'^2/c4'^2/ 
(v+l)/(k+l)$ 

As happened above, evaluating the resultant of t(l) and t{2) with respect to u is 
unworkable and we must use the same type of argument employed above in order to 
show that this resultant is not identically zero. Consider the particular value k — 2 
in t{l) and t{2). The resulting polynomials, called tt{l) and tt(2) arc shown to be of 
the same degree as t{l) and t(2) respectively (as polynomials in u). We then show 
that tt{l) and tt{2) have no common zeros. Hence the resultant of t{l) and t{2) with 
respect to u is not identically zero. 

tt(l):=sub(k=2,t(l))$ length coeff(tt(l),v); 9 length coeff(t(l),v); 9 
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tt(2):=sub(k=2,t(2))$ length coeff(tt(2),v); I4 length coeff(t(2),v); I4 
groebner{tt(l),tt(2)}; 

{1} 

Since the resultant(t(l),t(2),u) depends only on k and it is non-zero, we can con- 
clude that k must be constant. It only remains to show that v, u and s are also 
constants. This is shown in the next few lines 

groebner(coeff(t(l),v)); 

{1} So, V must be constant 
factorize lterm(h(l),u); 

{6 , u, u, u, k, k — 1 , k + 1 , 3k + 1} u must also be constant 
factorize Iterm(z3,s); 

{32, —s, s, k, k + 1 , 3k + 1} s is also constant. 

But k,s,u,v constants is impossible when gam(2,2,3)7^ (because, in particular, 
their time derivative should vanish and we have excluded this case before). Thus, 
it follows that gam(2,2,3) 7^ on is impossible. At the very beginning we also 
showed that gam(2,2,3) = and gam(l,l,3) 7^ on a non-empty open subset is im- 
possible. Hence, gam(l,l,3) = gam(2,2,3) = everywhere on W. Furthermore, there 
is an obvious symmetry between the coordinates x, y and z (there is nothing special 
about z). It follows necessarily gam(l,l,3) =gam(2,2,3) = gam(l,l,2) =gam(3,3,2) = 
gam(2,2,l) =gam(3,3,l) = everywhere on W. Since W is dense on Up x Iq (we have 
restricted W always to smaller open dense subsets thereof) all these Christoffel sym- 
bols vanish identically on Up Vt G Iq. In terms of the metric ( plf ) at t = 0, this means 
ds'^\t=o = A{x)dx'^ + B{y)dy^ + C{z)dz'^. Performing a trivial coordinate change we 
obtain ds'^\t=o = dx^ -\- dy"^ + dz^. The Codacci equations D^aK''^ = D^aS^^ = imply 
that i^^, are diagonal constant matrices in the coordinate system {x,y,z}. Then, 
the evolution equations (||), (0) imply that K^{t) and S^{t) are independent of 
z. Finally, the evolution equation (|^) shows that habit) is diagonal and depends only 
on t (in the coordinate system {x, 5}). After dropping the tildes we readily obtain 
ds'^{t) = A{t)dx^ + B{t)dy'^ + C{t)dz^, which was the aim of this Appendix. 
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